Predictions for the intrinsic UV continuum properties of 
star forming galaxies and the implications for inferring 
dust extinction 



Stephen M. Wilkins 1 *, Violeta Gonzalez-Perez 2 , Cedric G. Lacey 2 , Carlton M. Bau 

1 University of Oxford, Department of Physics, Denys Wilkinson Building, Keble Road, OX1 3RH, U.K. 

2 Institute for Computational Cosmology, Department of Physics, University of Durham, South Road, Durham, DH1 3LE, U.K. 



March 2013 



ABSTRACT 

The observed ultraviolet continuum (UVC) slope is potentially a powerful diagnostic of 
dust obscuration in star forming galaxies. However, the intrinsic slope is also sensitive 
to the form of the stellar initial mass function (IMF) and to the recent star forma- 
tion and metal enrichment histories of a galaxy. Using the GALFORM semi-analytical 
model of galaxy formation, we investigate the intrinsic distribution of UVC slopes. For 
star-forming galaxies, we find that the intrinsic distribution of UVC slopes at z = 0, 
parametcriscd by the power law index /3, has a standard deviation of ap ~ 0.30. This 
suggests an uncertainty on the inferred UV attenuation of ^4f uv ~ 0.7 (assuming a 
Calzetti attenuation curve) for an individual object, even with perfect photometry. 
Furthermore, we find that the intrinsic UVC slope correlates with star formation rate, 
intrinsic UV luminosity, stellar mass and redshift. These correlations have implica- 
tions for the interpretation of trends in the observed UVC slope with these quantities 
irrespective of the sample size or quality of the photometry. Our results suggest that 
in some cases the attenuation by dust has been incorrectly estimated. 
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1 INTRODUCTION 

The ultraviolet continuum (UVC) slope has been proposed 
as a robust diagnostic of dust attenuation in star forming 
galaxies (e.g. Meurer et al. 1999). The usefulness of the UVC 
slope is most apparent at high- redshift where alternative di- 
agnostics of dust attenuation such as the far-IR are generally 
inaccessible and galaxies are selected on the basis of their 
rest-frame UV emission (e.g. Wilkins et al. 2010; Wilkins et 
al. 2011a). This has prompted several studies which care- 
fully study the UVC slopes of high-redshift galaxies in an 
attempt to constrain the degree of dust obscuration in UV 
selected samples (e.g. Bouwens et al. 2009, Dunlop et al. 
2011, Wilkins et al. 2011b, Bouwens et al. 2011, Finkelstein 
et al. 2012). 

In all of these studies a simple linear relation is assumed 
between the observed UVC slope and the amount of dust at- 
tenuation. This requires the assumption of a unique intrinsic 
UVC slope. However, the intrinsic UVC slope is also affected 
by a number of other properties besides dust, including the 
recent star formation and metal enrichment histories, and 
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the form of the stellar initial mass function (IMF). The sen- 
sitivity of the intrinsic slope to these properties potentially 
limits the usefulness of the observed UVC slope as an ac- 
curate diagnostic of dust attenuation. In this work we em- 
ploy the GALFORM semi-analytical model of galaxy forma- 
tion (first developed by Cole et al. 2000) to produce realistic 
star formation and metal enrichment histories with which to 
determine both the intrinsic distribution of UVC slopes j3i 
and to investigate how this distribution varies with stellar 
mass, star formation rate, UV luminosity and redshift. Lacey 
et al. (2011) have shown this model can reproduce the ob- 
served UV luminosity function over a wide range of redshifts 
(see also Gonzalez- Perez et al. 2012, in prep). 

This paper is organised as follows. In i|2.1l and i|2.2l 
we introduce the GALFORM semi-analytical galaxy forma- 
tion model and filter convention respectively. In Section [3] we 
discuss the principal physical properties which affect the ob- 
served UVC and investigate, using the PEGASE.2 stellar pop- 
ulation synthesis (SPS) model (Fioc & Rocca-Volmerange 
1997, Fioc & Rocca-Volmerange 1999), how the UVC slope 
is affected by dust and changes in the recent star formation 
and metal enrichment histories, and by the choice of IMF. 
In Section [3] we also investigate the distribution of UVC 
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colours/slopes predicted by the Galform model as a result 
of the predicted variety in the star formation history (SFH) 
alone (fj3T2j) and in the metal enrichment history (i )3.3p . In 
Section [4] we investigate how the distribution of intrinsic 
UVC slopes is affected by the star formation rate f i]4.1[l . in- 
trinsic UV luminosity ( H4.2p , stellar mass ( H4.3P and redshift 
( i]4,4|l . In Section[5]we present our conclusions. 



2 MODELLING APPROACH 

In this section we give an overview of the GALFORM galaxy 
formation model f i)2.1[) and set out how the UV continuum 
slope is measured ( H2.2I 1. 

2.1 The Galaxy Formation Model 

We predict the intrinsic UV-slope of galaxies in a ACDM 
universe using the GALFORM semi-analytical galaxy forma- 
tion model developed by Cole et al. (2000). Semi-analytical 
models use physically motivated equations to follow the fate 
of baryons in a universe in which structure grows hierarchi- 
cally through gravitational instability (see Baugh 2006 for 
an overview of hierarchical galaxy formation models). 

GALFORM follows the main processes which shape the 
formation and evolution of galaxies. These include: (i) the 
collapse and merging of dark matter haloes; (ii) the shock- 
heating and radiative cooling of gas inside dark matter 
haloes, leading to the formation of galaxy discs; (iii) qui- 
escent star formation in galaxy discs; (iv) feedback from su- 
pernovae, from active galactic nuclei (AGN) and from pho- 
toionization of the intergalactic medium (IGM); (v) chemi- 
cal enrichment of stars and gas; (vi) galaxy mergers driven 
by dynamical friction within common dark matter haloes, 
leading to the formation of stellar spheroids, which also may 
trigger bursts of star formation. The end product of the cal- 
culation is a prediction for the number and properties of 
galaxies which reside within dark matter haloes of different 
masses. The outputs include the star formation and metal 
enrichment histories for each galaxy, including the contri- 
bution of merger driven starbursts. The star formation his- 
tories predicted by the model are in general more complex 
than the exponentially decaying models typically assumed 
to interpret observations (see Baugh 2006 for examples). 

In this paper we focus our attention on the Baugh et al. 
2005 (hereafter B05) model. Some of the key features of this 
model are (i) a time-scale for quiescent star formation that 
is assumed to vary as a power of the disc circular velocity 
(see Lagos et al. 2011 for a study of different star forma- 
tion laws in quiescent galaxies) , (ii) bursts of star formation 
are triggered only by galaxy mergers, (iii) the default im- 
plementation of this model adopts a Kennicutt (1983) IMF 
(£ = dN/dm oc m~ 25 for m > 1M Q and £ oc m~ 1A for 
m < IMq, c.f. Salpeter 1955: £ oc m~ 2,35 ) in quiescent star 
formation in galactic discs, while in starbursts a top-heavy 
IMF (£ oc m -1 ) is assumed, (iv) the inclusion of super- 
nova feedback with superwinds (see Benson et al. 2003 for 
a discussion of the effect that feedback has on the lumi- 
nosity function of galaxies), and (v) the reionisation of the 
intergalactic medium is approximated by a simple model in 
which gas cooling is completely suppressed in haloes with 
circular velocities less than 30kms _1 at redshifts z < 10 



(Lacey et al. 2011). The parameters of this model were 
fixed with reference to a subset of the available observations 
of galaxies, mostly at low redshift. The B05 model uses, 
by default, the simple stellar population (SSP) spectral en- 
ergy distributions (SEDs) generated by Bressan, Granato, & 
Silva (1998), using the Padova 1994 stellar evolution tracks 
and the model stellar atmospheres from Kurucz (1993)0- The 
B05 model uses the canonical (ACDM) parameters: matter 
density, flo = 0.3, cosmological constant, — 0.7, baryon 
density, ilb = 0.04, a normalisation of density fluctuations 
given by ag = 0.93 and a Hubble constant today Ho = 70 
kms _1 Mpc _1 . The B05 model employs merger trees gener- 
ated using the Monte Carlo algorithm introduced by Parkin- 
son et al. (2008). The B05 model also includes a dust model 
in which stars and dust are mixed together. The dust ex- 
tinction in GALFORM and the impact on the rest-frame UV 
is discussed in Lacey et al. (2011) and Gonzalez-Perez et al. 
(2012, in prep). In the calculations presented here we study 
the intrinsic UVC predicted by GALFORM, and so omit dust 
extinction in the model predictions. 

In addition to reproducing local galaxy data, the B05 
model matches the number and redshift distribution of 
galaxies detected by their emission at sub-millimetre wave- 
lengths, the luminosity function of Lyman break galaxies 
and the abundance and clustering of Lyman-alpha emitters 
(Orsi et al. 2008). No parameters have been tuned for the 
study presented here. We refer the reader to B05 and Lacey 
et al. (2008, 2011) for a full description of this model. As a 
method of exploring differences in the distribution of UVC 
slopes arising from the initial mass function, star formation 
history and metal enrichment history we consider 3 seperate 
implementations of the B05 model: 

• default - the default implementation of the model as 
described by B05. 

• single IMF - same as the default implementation except 
the top-heavy IMF in starbursts is replaced by the Kennicutt 
IMF, so that this IMF is used in all modes of star formation. 

• single IMF + single metallicity - same as the default 
implementation except the top-heavy IMF in starbursts is 
replaced by the Kennicutt IMF and all stars are assumed to 
have the same metallicity (Z = Zq). 



2.2 Filters 

In order to reproduce the techniques commonly applied to 
observations, we measure the UVC slope using two artificial 
broadband rest-frame filters (shown in Fig. [TJ, FUV {T\ — 
[0.13 < X/fim < 0.170.) and NUV (T\ = [0.18 < X/fim < 
0.26]). These filters cover a similar wavelength range to the 
GALEX filter set. The use of rest-frame filters allows us 
to consistently compare the distribution of UVC slopes in 
galaxies at different redshifts. 



1 These are the same isochorones and stellar atmospheres used 
by Bruzual & Chariot (2003). 

2 We utilise the Iverson bracket notation such that [A] = 1 when 
A is true and otherwise. 
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Figure 1. The synthetic UV spectrum of a star forming galaxy 
(100 Myr of previous continuous star formation, Salpeter IMF, 
Z = Zq) produced with the PEGASE.2 population synthesis model 
(black line). The heavy grey line is a pure power law fx oc A^ as- 
suming the value of /? inferred from the FUV-NUV colour (trans- 
mission functions shown in upper panel). 

2.2.1 Parameterisation of the UVC slope 

Observations of the UVC colours of galaxies at different red- 
shifts are difficult to compare as in each case the available 
filters probe different ranges of the UVC. To compare the 
UVC properties at different redshifts it is instead common 
to parameterise the continuum by a power law, with slope 
P i.e.: 

/aocA", (1) 
or in terms of /„, 

/, cc A 2 +^ . (2) 

In the context of our choice of rest-frame filters, the conver- 
sion from (FUV — NUV)ab colour to /3 is: 

P = C x (FUV - NUV)ab - 2 , (3) 

where C is a conversion factor determined by convolving the 
two filter transmission functions with a power law spectrum. 
The value of C is sensitive to the individual filters and for 
the FUV and NUV filters defined above is C = 2.49. Fig. 
[1] shows both the synthetic UV spectrum of a star forming 
galaxy and the power law using the value of /3 inferred from 
the FUV - NUV colour. 

Throughout this work we make a distinction between 
the intrinsic UVC slope /3i and the observed slope f3 . The 
intrinsic slope is assumed to be the slope produced by stars 
in the absence of any extrinsic effects such as dust extinction 
or nebular emission. 



3 PHYSICAL PROPERTIES AFFECTING THE 
OBSERVED UV CONTINUUM SLOPE 

The observed ultraviolet continuum slope of a galaxy is af- 
fected by various properties both intrinsic and extrinsic to 
the stars producing the UV emission. 

3 For a combination of observer frame niters C is also sensitive 
to the source redshift. 



The intrinsic UVC of actively star forming galaxies 
is dominated by emission from the most massive, high- 
temperature stars (m > 2 Mq , OB A class) . The UV 
(1216 — s> 3000A) spectral flux density of these stars is 
well characterised by a power law with a blue slope (i.e. 
if /„ « \ 2+f) then < -2) as shown in Fig. [TJ This power- 
law behaviour principally arises because the UV continuum 
probes the Rayleigh- Jeans tail of the black body spectral en- 
ergy distribution (though it deviates from a pure Rayleigh- 
Jeans behaviour because of the effects of opacity in the star's 
atmosphere) as shown by the synthetic spectrum in Fig. 
[T] At lower temperatures, the range 1216 — > 3000A probes 
where the black body distribution peaks. This effect, com- 
bined with variation in the opacity (as a function of tem- 
perature etc.), results in the SEDs of the cooler, lower-mass 
stars having a redder UV continuum than those of stars with 
higher temperatures/masses. The sensitivity of the UVC of 
a star to its temperature (and opacity), and thus metallic- 
ity, mass and age means that the intrinsic continuum of a 
composite stellar population is sensitive to the distribution 
of stellar masses, ages and metallicities. The distribution of 
masses is in turn determined by the star formation history 
and IMF; thus both these factors affect the intrinsic UV 
continuum colours. 

In addition to these intrinsic properties the observed 
UVC is also affected by the presence of intervening dust 
through scattering and absorption and (typically to a much 
lesser extent) continuum emission from hot ionised gas. 

In the following sections (§3.1-3.4) we both outline the 
effect of dust (§3.1) and investigate the effect of intrinsic 
properties (including the star formation history, metallicity 
and IMF) on the observed UVC colours/slope. 



3.1 Dust 

The principal property which causes the observed UVC slope 
Po of luminous star forming galaxies to differ from the intrin- 
sic slope is dust. Both empirically determined (e.g. Cardelli 
et al. 1989, Calzetti et al. 2000 - hereafter COC0), and the- 
oretically motivated dust attenuation laws suggest that at- 
tenuation increases rapidly from the NIR/optical through 
the UV. As such, the UV colour is potentially a powerful 
diagnostic of dust attenuation as the presence of dust will 
redden the observed slope compared to its intrinsic value. 
The exact relation depends on the characteristics of the at- 
tenuation curve in the UV; for example, depending on the 
choice of filters and redshift, a curve including the bump 
feature at A = 2175A (such as the Cardelli et al. 1989 ex- 
tinction curve) implies a much weaker relationship between 
the UVC slope and dust than the COO attenuation curve 
which decreases monotonically with wavelength. 

Here we use the COO attenuation curve to see the effect 
of dust on the observed UVC slope. Assuming the underlying 
intrinsic UVC is described by a power law (Eq.[5J) a relation- 
ship between the observed and intrinsic UVC slopes, p o and 



4 The Calzetti et al. (2000) starburst attenuation curve encapsu- 
lates the effects of scattering and absorption while the Cardelli et 
al. (1989) curve describes extinction by a foreground screen and 
includes the prominent bump at A = 2175A. 



4 Stephen M. Wilkins, et al. 



Pi respectively, and the attenuation in the FUV, Af nv , can 
be derived (c.f. COO), 



A fuv = 2.37 x [/3o - Pi]. 



(4) 



Assuming a synthetic intrinsic UVC (produced in this ex- 
ample using the PEGASE.2 SPS model) corresponding to 
100 Myr previous constant star formation, and a Salpeter 
IMF with solar metallicity (as shown in Fig. [l]) yields an 
intrinsic UVC slope of Pi ~ —2.3. Inserting this into Eq. [4] 
yields 



A fuv = 2.37/3 + 5.5, 



(5) 



which is similar to the relation proposed by Meurer et al. 
(1999). This can be rewritten to express the uncertainty on 
the attenuation arising from the uncertainty in the intrinsic 
UVC slope, 



5A fuv = 2.37 x 



(6) 



This highlights the strong dependence of the inferred UV 
attenuation on the intrinsic UVC slope. 
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Figure 2. The effect of the duration and form of previous star 
formation on the UVC colour. The three lines show the evolution 
of the UVC colour (left) and slope (right) of a stellar population 
forming stars at a constant rate (solid line), an exponentially in- 
creasing rate (dashed line) and as an instantaneous burst (dotted 
line). 



3.2 Star Formation History 

The strong variation with mass of the main-sequence life- 
times of stars means that the contemporary mass distribu- 
tion, and thus UV continuum slope, is sensitive to a galaxy's 
star formation history. After less than 10 Myr of continuous 
star formation, a stellar population would still retain most of 
its original high-mass stars (i.e. up to this time, the contem- 
porary mass distribution is similar to the IMF). However, af- 
ter prolonged periods of star formation some fraction of the 
most massive stars will have evolved off the main-sequence. 
This reduces the relative contribution of these stars to the 
UV continuum, resulting in a redder slope (i.e P > —2). 

The limiting cases of the effect of the previous star for- 
mation history are shown in Fig. [2] Here the evolution of 
the UVC slope is shown for three cases, an instantaneous 
burst, a constant SFH and an exponentially increasing SFH 
(assuming solar metallicity and a Salpeter IMF). The last 
case essentially leaves the slope constant, as the UVC is con- 
tinuously dominated by the most massive OB-type stars. In 
the first case, the intrinsic UVC slope quickly reddens as it 
becomes dominated by progressively cooler stars. For the in- 
termediate case, corresponding to a constant star formation 
rate, the intrinsic UVC slope evolves mildly, changing by 
only 8 Pi — 0.3 between star formation durations of 10 Myr 
and 1 Gyr. 

Fig. [2] highlights that the UVC slope is sensitive to the 
recent star formation history. Estimating the uncertainty on 
the intrinsic UVC slope, 6 Pi , requires knowledge of the range 
of star formation histories. To obtain these we use simulated 
galaxies taken from the GALFORM model. In order to inves- 
tigate the effect of the star formation history on the UVC 
we consider special cases in which the model is run with a 
fixed (solar) metallicity and a single IMF (Kennicutt 1983) , 
with the variation in the intrinsic UVC slope then being en- 
tirely due the effect of the star formation history. The result- 
ing distribution of UVC slopes, for galaxies selected as star 
forming (i.e. with star formation rates ip > O.lM0yr _1 ), is 
shown in Fig. [3] The structure of this distribution is asym- 
metric with a long red tail. The asymmetry in the distri- 
bution (and the long red tail) is driven by the fact that 
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Figure 3. The simulated distribution of intrinsic UVC slopes (top 
axis) and colours (bottom axis) for star forming galaxies (defined 
as those with SFR> 0.1 Moyr -1 ) in the GALFORM model at z = 0, 
assuming a single constant metallicity and universal IMF. The 
solid, dotted and dashed vertical lines denote the weighted mean, 
median and the 15.87th to 84.13th percentile range respectively. 
The solid horizontal line shows the standard deviation around the 
median. This, and subsequent, histograms are normalised by the 
modal value. 



while it is possible to obtain extremely red UVC slopes by 
having a dominant old stellar population, it is impossible 
to push the UVC slope to be arbitrarily blue. Despite the 
non-gaussian structure of this distribution, the standard de- 
viation around the median provides a close indication of the 
15.87th to 84.13th percentile confidence interval. The stan- 
dard deviation of this distribution is ap ~ 0.13; assuming 
the formula relating the uncertainty on Pi to the UV at- 
tenuation (Eq. [6}, this suggests an intrinsic uncertainty on 
the UV attenuation due to the variation in the SFH of the 
galaxy population of 5A{ UV ~ 0.31. 
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Figure 4. The effect of assuming different values of the stellar 
metallicity on the intrinsic UVC colour (left) and slope (right) 
derived from the Pegase.2 SPS model (assuming assuming a 
100 Myr previous duration of continuous star formation prior to 
measurement and a Salpetcr IMF). 



3.3 Metallicity 

Both the opacities and effective temperatures, and thus the 
ultraviolet continua, of stars are affected by their chemical 
composition; stars with lower metal abundances typically 
have bluer UVC slopes. Hence, the intrinsic UVC slope is 
also sensitive to the metallicity of the UV luminous pop- 
ulation, or specifically to the recent metal enrichment his- 
tory. Fig. [4] shows the impact of metallicity on the intrinsic 
UVC slope (assuming constant star formation in the 100 Myr 
prior to observation and a Salpeter (1955) IMF, using the 
PEGASE.2 population synthesis model); changing the metal- 
licity from Z — 0.02 to 0.004 causes the intrinsic UVC slope 
Pi to steepen by ~ 0.25. Thus the UV attenuation, Af uv , 
of a galaxy which has Z — 0.004, but which is assumed to 
have Z = 0.02 would be overestimated by ~ 0.6 mag due to 
the difference in the intrinsic UVC colour. 

As with the SFH, accounting for the metallicity effects 
is difficult without some knowledge of the metal enrichment 
history. Once again we turn to the GALFORM model to pro- 
vide plausible distributions of star formation and metal en- 
richment histories from which we can determine the intrinsic 
UVC slope. Including the predicted variation in metal en- 
richment histories of galaxies widens the intrinsic distribu- 
tion of UVC slopes compared to that obtained assuming a 
single metallicity (the resulting distribution is shown in Fig. 
[5]l . The standard deviation of the intrinsic UVC slope dis- 
tribution increases to crp ~ 0.30 compared with op ~ 0.13 
with no metallicity variation. Using Eq. [6l this suggests an 
intrinsic uncertainty in the UV attenuation inferred from the 
UVC colour of 8Af uv ~ 0.7. This suggests that, even in the 
absence of photometric noise or redshift uncertainties, the 
dust attenuation of a single object can not be measured more 
accurately than 5Af uv ~ 0.7 without additional information 
to constrain the star formation history or metallicity. 



3.4 Initial Mass Function 

The IMF (£(m) = AN /Am) describes the stellar mass distri- 
bution of a single stellar population with £ agc = 0. A range 
of parameterisations of the IMF exist, with the simplest tak- 
ing the form of single power law (i.e. £(m) oc m a ), e.g. the 
Salpeter (1955) IMF with a — —2.35. Updated parameteri- 
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Figure 5. The simulated distribution of UVC colours for star 
forming galaxies (SFR> 0.1 M0yr _1 ) at z = assuming a single 
initial mass function but varying metal enrichment and star for- 
mation history. The solid, dotted and dashed vertical lines denote 
the weighted mean, median and 15.87th and 84.13th percentile 
range respectively. The solid horizontal line shows the standard 
deviation around the median. 



0.2 



-0.2 



-0.4 



-2 

-2.2 

-2.4 

-2.6 

-2.8 



-1.5 -2 

a 2 (high-moss slope) 



-2.5 



Figure 6. The effect of the choice of the high-mass slope «2 of 
the IMF on the UVC colour/slope derived from the Pbgase.2 
SPS mode (assuming assuming a 100 Myr previous duration of 
continuous star formation and solar metallicity). 



sations take into account the observed flattening of the IMF 
below some characteristic mass (e.g. m c ~ 1 Mq for the Ken- 
nicutt IMF) by adopting a broken power law parameterised 
by ai and cii (the low and high-mass slopes respectively). 

The IMF can be modified in a number of ways in the 
context of this parameterisation. One method is to simply 
vary the high-mass slope cti. Steepening the high-mass slope 
(decreasing a.2) will reduce the relative contribution of very- 
high mass stars to the UV luminosity; this will have the ef- 
fect of reddening the UVC colour/slope. This can be seen 
in Fig. [6] (assuming 100 Myr previous duration of continu- 
ous star formation and solar metallicity); flattening «2 from 
—2.35 — > —1.5 reduces the slope by A/3j ~ 0.2 assuming the 
same star formation and metal enrichment history. 

The effect of the IMF is potentially important as a num- 
ber of studies have suggested the high-mass slope of the IMF 
may be flatter than Salpeter (e.g. Wilkins et al. 2008b) or 
that the IMF may effectively vary from galaxy to galaxy (e.g. 
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Baldry & Glazebrook 2003, Wilkins et al. 2008a, Hoversten 
& Glazebrook 2008, Lee et al. 2009, Finkelstein et al. 2011). 

The default implementation of the B05 GALFORM model 
adopts a top-heavy IMF (£ oc m _1 ) in merger triggered SF, 
while retaining the Kennicutt IMF in quiescent SF. At low 
redshift, where our attention has been focused thus far, the 
fraction of star formation occurring with the flat, top-heavy 
IMF is small in the B05 model. This fraction increases with 
redshift such that at z ~ 3 the contribution to the star for- 
mation rate density from each mode is roughly similar. In 
Fig. we compare the intrinsic UVC slope distribution of 
galaxies at z ~ 3 using both the default implementation of 
B05 model (i.e. with the top-heavy IMF mode of star for- 
mation in bursts) and a variant in which a universal IMF is 
used. The main difference between the two is that the distri- 
bution of UVC colours in the default B05 model is on average 
bluer and the width of the distribution slightly smaller. A 
second peak also emerges associated with burst-driven SF. 
This can be seen clearly in Fig. [7] where the distribution 
for galaxies whose total star formation is dominated by the 
burst mode is shown alongside the full distribution. 



4 TRENDS WITH SFR, UV LUMINOSITY, 
STELLAR MASS AND REDSHIFT 

In the previous section we showed that the distribution of 
intrinsic UVC slopes (of star forming galaxies) is affected by 
the IMF, and the star formation and metal enrichment histo- 
ries of galaxies. In this section we investigate if the median 
intrinsic UVC slope correlates with a galaxy's star forma- 
tion rate, intrinsic UV luminosity, stellar mass and redshift. 
A correlation of the intrinsic slope with any of these proper- 
ties may result in an erroneous observed correlation between 
the property and the inferred dust extinction, or may rein- 
force a weak or negligible correlation. 

Fig. [8] shows the median intrinsic UVC slope for the 
three variants of the B05 model (single IMF and metallicity, 
single IMF, and default) binned by the star formation rate, 
intrinsic UV luminosity and stellar mass, which we discuss 
in turn below. 
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Figure 7. The simulated distribution of intrinsic UVC colours 
and slopes for star forming galaxies (tp > 0.1 MQyr - 1 ) at z = 3 
assuming the B05 model with a single Kennicutt (1983) IMF 
(top) and the default implementation with two IMFs (bottom). 
The histogram enclosed by the dark line shows the distribution 
for galaxies whose total SF is dominated by the burst mode while 
the hatched histogram shows those dominated by quiescent star 
formation (the grey shaded histogram is for all star forming galax- 
ies). 



4.1 Star Formation Rate 

The top panels of Fig. [S] show the median UVC slope (and 
15.87th to 84.13th percentile range) for galaxies binned by 
their star formation rate at z = and z = 6 using the three 
variants of the B05 model (single-metallicity and single- 
IMF, single-IMF and default). For the single-metallicity sin- 
gle IMF model there is no correlation between the star 
formation rate and median intrinsic UVC slope, irrespec- 
tive of redshift. However for both the single-IMF and de- 
fault variants of the model there is a negative correlation 
at low- redshift. This trend becomes progressively flattened 
at higher-redshift and for z > 4 is essentially flat (z = 6 is 
shown in Fig. [HJ. The existence of this trend in the single- 
IMF and default variants of the model and its absence in the 
single-metallicity and single-IMF model suggests this trend 
is driven by the metallicity variation as a function of star 
formation rate. At high-redshift the variation in metallicity 
is less pronounced leading to the absence of any correlation. 



4.2 Intrinsic UV Luminosity 

Related to the star formation rate but more useful in an 
observational context is the intrinsic FUV luminosity, Lf uv . 
The middle panels of Fig. [5]shows the median intrinsic UVC 
slope binned by the intrinsic absolute FUV magnitude. Due 
to the close connection between the SFR and the intrin- 
sic UV luminosity a similar trend exists to that described 
above; at low-redshift there is a mild trend such that more 
luminous galaxies have an intrinsically bluer UVC colour (in 
the single-IMF and default variants of the model), while at 
higher-redshift the correlation is virtually flat. The fiat trend 
at high-redshift, reinforces the conclusions of several studies 
(e.g. Wilkins et al. 2011b) that galaxies with brighter ob- 
served UV luminosities have greater dust attenuation than 
their low luminosity counterparts. 
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Figure 8. The median intrinsic UVC slope binned by star formation rate (top), intrinsic FUV absolute magnitude (middle) and stellar 
mass (bottom) at z = (black) and 2 = 6 (grey) for each of the three model variants: single-metallicity and single-IMF (left), single-IMF 
(middle) and default (right). In each case the horizontal bar denotes bin width while the vertical bar denotes the 15.87th to 84.13th 
percentile range of the GALFORM model predictions. 



4.3 Stellar Mass 

The lower panel in Fig. [8] shows the variation of ft with stel- 
lar masses for galaxies around M* at z — and z — 6. At 
z — (but extending to 2 ~ 1), and for all three model vari- 
ants, the median UVC colour is also correlated with stellar 
mass, becoming redder at high-stellar masses. The width (as 
measured by the confidence interval) of the UVC slope dis- 
tribution also increases dramatically to high stellar masses. 
The fact that this trend exists for all three models suggests 
that it is driven predominantly by the star formation history 
in so far as more massive galaxies are typically composed of 



older populations (e.g. Kauffmann et al. 2003). At higher- 
redshift this correlation flattens, to the extent that at z — 6, 
the trend is virtually flat. 

One implication of the strong correlation between stel- 
lar mass and UVC slope is that assuming a constant intrinsic 
slope to convert the observed slope to an attenuation may 
introduce an artificial correlation. As an example, apply- 
ing the same intrinsic slope found at log 10 (Mt/M©) = 9 
(Pi ~ —2.2) would result in a systematic over-estimation 
of the attenuation at log 10 (M* /Mq) — 11 (/3j ~ —1.6) of 
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Figure 9. The evolution of the median UVC slope of star forming 
galaxies as a function of redshift for both the single IMF and 
default implementation of the B05 model. 

SAf uv ~ 1.42 assuming the single-IMF variant of the B05 
model. 

4.4 Redshift 

Potentially the most important trend is the evolution of the 
median UVC slope with redshift. Shown in Fig. [9] is the 
redshift evolution of the median UVC slope for star forming 
galaxies (SFR> 0.1 MQ-yr^ 1 ) assuming both the single IMF 
(solid line) and default (dashed line) implementations of the 
B05 model. 

In both cases the median intrinsic UVC slope evolves, 
becoming progressively redder towards lower redshift. The 
cause of this is a combination of both the increasing average 
metallicity of recently formed stars and the changing stellar 
mass distribution in star forming galaxies. As progressive 
generations of stars form, the mass distribution in galaxies 
becomes increasingly dominated by low and intermediate 
mass stars with long main-sequence lifetimes. The evolution 
of the median UVC slope assuming the default implementa- 
tion of the B05 model is similar to the single-IMF implemen- 
tation of the B05 model, though slightly bluer over the entire 
history. Specifically, at z > 1 the difference between models 
is significant (8 Pi — 0.1) while at lower-redshift (z < 0.5) 
the difference becomes negligible 8 Pi < 0.05, reflecting the 
smaller contribution of merger-driven star formation to the 
total star formation rate density at these epochs. 

If this decrease in the median intrinsic UVC slope is not 
properly taken into account, and, for example, a relation rel- 
evant at low-redshift is applied to high-redshift galaxies (as 
is typically done in the literature), the UV attenuation in- 
ferred will be systematically underestimated. For example, 
applying a relation based on the intrinsic UVC slope dis- 
tribution at z = to galaxies at z = 7 would result in 
the systematic under-estimation of the UV attenuation by 
8Af uv ~ 0.97 and the possible introduction of an erroneous 



redshift trend or at least an over-estimation of the redshift 
evolution. 

In addition, in both cases, the width of the distribution 
also evolves, becoming slightly narrower at higher redshift 
(at z = 6 op = 0.26, c.f. z = 6 op = 0.30 for the single-IMF 
implementation) . For the default implementation we also see 
the emergence of bimodality in the distribution as a result of 
the merger-driven SF occurring with a top-heavy IMF. At 
low redshift, where the contribution of merger-driven star 
formation is small, and at very-high redshift, where the con- 
tribution of quiescent star formation is low, the bimodality 
essentially disappears. 



5 CONCLUSIONS 

We have used realistic star formation and metal enrichment 
histories predicted by the GALFORM semi-analytical galaxy 
formation model to investigate the intrinsic distribution of 
ultraviolet continuum (UVC) colours and slopes. We find, 
for star forming galaxies at low-redshift, that the standard 
deviation of the distribution of UVC slopes (parameterised 
by the power-law index Pi) is op ~ 0.30 (assuming a single 
IMF). Assuming the Calzetti et al. (2000) reddening curve 
this suggests an intrinsic uncertainty of 8A1500 ~ 0.71 in the 
UV attenuation inferred from the UVC slope for an individ- 
ual object without any aditional information to constrain 
the star formation or metal enrichment histories. 

We also investigated how the median intrinsic UVC 
slope and the width of the pi distribution correlates with 
various properties including the star formation rate, intrinsic 
UV luminosity, stellar mass and redshift. At low-redshift we 
find that the median UVC slope is sensitive to the star for- 
mation rate, intrinsic UV luminosity and stellar mass, with 
more massive galaxies typically having redder UV contin- 
uum slopes and those with higher star formation rates (and 
higher UV luminosities) typically being bluer. At higher- 
redshift these correlations flatten and become less impor- 
tant. The model also suggests significant evolution of me- 
dian UVC slope with redshift (though the width of the dis- 
tribution changes only mildly). These various correlations 
suggest that trends of the observed UVC slope with stellar 
mass or redshift are not driven entirely (if at all) by dust 
but partly by the evolution of the star formation and metal 
enrichment histories of galaxies. In order to correctly inter- 
pret observations, in the context of dust reddening, it is then 
important to compare the observed distribution of UV con- 
tinuum slopes with a model including a physically motivated 
treatment of dust but taking also into account the predicted 
intrinsic distribution. In a follow up paper (Wilkins et al. 
2012) we compare the observed distribution with the pre- 
dicted intrinsic one for bright z ~ 3 — 5 star forming galax- 
ies. 
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